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Abstract 

From a fit to the experimental data on the 66 fine structm'e, the two-loop strong 
coupling constant is extracted. For the IP state the fitted value is asi/J-i) = 
0.33 ± 0.01 (exp) lb 0.02 (th) at the scale /ii = 1.8 ± 0.1 GeV, which corresponds to the 
QCD constant A^"^) (2 - loop) = 338 ± 30 MeV (n/ = 4) and as{Mz) = 0.119 ± 0.002. 
For the 2P state the value as(/i2) = 0.40 ± 0.02 (exp) ± 0.02 (th) at the scale 
1^2 = 1.02 lb 0.02 GeV is extracted, which is essentially larger than in the previous 
analysis of refs. |Q, ^, but about 30% smaller than the value given by the standard 
perturbation theory. This value 0^(1.0) ~ 0.40 can be obtained in the framework of 
the background perturbation theory thus demonstrating the freezing of as(/i). The 
relativistic corrections to are found to be about 15%. 

PACS 11.10.Ij,11.15Bt,12.38.Lg 
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1 Introduction 



The bottomonium spectrum is one of the richest among all known mesons and its levels were 
measured with high precision [|l|. These data about bb states have been intensively studied 
in different theoretical approaches, in particular, to determine the QCD strong coupling 
constant at different energy scales /i from the level differences [2-10]. At present, 

however, there is no clear picture which are the exact values of as{fi) for the bb levels and 
how they are changing from the ground state to the excited ones. There are several reasons 
for that. 

First of all, there is no experimental information on the r]b{nS) masses and therefore 
Osifi) cannot be directly determined from the bb hyperfine splittings in S-wave states. 

Second, to describe the fine structure splittings in the P-wave states, different energy 
scales were used in different theoretical analyses [4-7]. In ref. [Q ^^(yu) = 0.33 (yU = 
3.25 GeV) was taken for all bb S- and P-wave states, while in [^ fi was chosen to be equal 
to the b quark mass, fi = m with either m = 4.6 GeV or m = 5.2 GeV. The fitted values of 
as{fi) were found to be as{m) = 0.22 -i- 0.27 [^ and for the 2P state appeared to be 

smaller than for the IP state. 

An important step to clarify this problem was taken in |^ where the low-lying bot- 
tomonium states, IS, 2S, and IP were investigated. It was observed there that the scale fj, is 
a decreasing function of the principal quantum number n,fi = 2{na)~^ where a is a Coulomb 
type radius. Therefore, /i is found to be equal for the 25* and IP states and the values 
fi = 1.7GeV, as(1.7) = 0.29, were determined from the fine structure splittings of Xfe(l-P)- 
Also, as(yu) is larger for excited states with a larger radius of the system, thus indicating that 
for a bound state the characteristic scale /i is determined by the size, but not by the mo- 
mentum of the system. One of our main goals here is to check this important statement for 
the 2P state, Xb(2P), which cannot be studied in the framework of the approach developed 
in 1,0. 

In the present study of the IP and 2P bb states we shall try to answer the following 
questions: 

What are the values of for the 2P and the IP states? 

Do the extracted values of correspond to the existing experimental data 

on asiMz) and A^"/)? 

How does as{fi) depend on the relativistic corrections to the wave functions in 
bottomonium? 

How sensitive are the extracted values of as{fi) to the b quark pole mass and the 
parameters of the static interaction? 

2 Perturbative Radiative Corrections 

It is well known that one cannot describe the fine structure splittings in heavy quarkonia 
without taking into account the second order radiative corrections [4-7, 10]. In coordinate 
space, perturbative static and spin-dependent potentials in the MS renormalization scheme 
were obtained in [^ ^ . From the potentials given there one can immediately find the matrix 
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elements of the spin-orbit and the tensor potentials: a = {VLs{r)), c = {Vrir)). Below we 
give their expressions for a number of flavours = 4, valid for the bb system: 
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and for the perturbative part of the the tensor splitting cp, 

Cp 



Cp -|- Cp , 



.(1) 



4 



.(2) 



4 ag(/i) 
3 7rm2 



/25 



7 



(r — In— + 5 +-(r In mr) 



V 6 



m 



6 



(1) 

(2) 

(3) 

(4) 



Here the constant A = f-fE + i = 1-44508 and = + f| = 3.42342. 

For our analysis it is convenient to introduce a linear combination of the matrix elements 

V 



a and c as was done in ref. 10 



|c — a. Its perturbative part r]p is 



n -^r n - ^r^^) ^(2) _ 
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The factor /4 in Eq. (5) can be found from Eqs. (2) and (4), 



h{nP) = - [1.97834 (r-=')„p - {r~Hnmr)nP 



(5) 
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For the fine structure analysis it turns out to be very important that the combination of 
matrix elements does not depend on the energy scale fi. Later, it will be also shown that 
/4 has the largest relativistic correction (about 35%) compared to other matrix elements and 
depends weakly on the parameters of the static interaction and on the mass of the 6-quark. 



3 Nonperturbative Contributions 

Besides the perturbative terms Eqs. (0, ^ the tensor and spin-orbit splittings have in general 
nonperturbative contributions: a = ap + Onp and c = Cp -|- Cnp. 

The nonperturbative part of the spin-orbit parameter onp is determined by the chosen 
confining potential. Here we shall take a linear potential ar at all distances. The linear 
behaviour of the confining potential at large distances is well established phenomenologi- 
cally due to the existence of the meson Regge trajectories and was also deduced from the 



minimal area law for Wilson loops in QCD (see the review [|TT[] and the references therein). 
Recently, the linear behaviour of the confining potential was found at small distances due 
to the interferences of perturbative and nonperturbative (NP) effects [|1^ and to the satu- 
ration property of the QCD strong coupling constant in vacuum field [13-15]. For the linear 
potential ar the nonperturbative interaction is given by the Thomas potential for which 

"NP = -^(r-> (7) 
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The nonperturbative contribution to the tensor sphtting can be found from the vacuum field 
correlator Di{x) |16[ which was measured in lattice QCD and was found to be of 



exponential form. Then, as was shown in |10|, |18|, [19 



c^P = §^^{r'Mr/T,))^^J{T,), J{T,) ^ ^^{r'Ko{r/T,)), (8) 

where Tg is the vacuum correlation length. Lattice QCD without dynamical fermions give 
Tg ^ 0.2 fm and Tg ^ 0.3 fm in the presence of dynamical fermions with four flavours [p!7[]. 

In refs. the correlator -Di(O) in Eq. (|^) was shown to be small: lattice calculations in 
quenched SU(3) theory give Di{0)/ D{0) ^ | and in full QCD with four staggered fermions 
i^i(0)/Z^(0) ~ 0.1, where D{0) is another vacuum field correlator which mostly determines 
the confining potential. These two correlators at the point x = can be expressed through 
the vacuum gluonic condensate G2 (here the vacuum correlators are normalized as in [p!l|, p!6| ) : 

D(0) + Di(0) = ^G2. (9) 

Therefore, the lattice estimate for Di{0)/D{0) is 0.1 -^ 0.3 and from the relation (9) one 
obtains 



^^^G2 < D,{0) < - G2. (10) 

Our calculations give the following typical values for the matrix elements J{Tg): for the IP 
state J{Tg) ^ 0.17 GeV"^ and for the 2P state J{Tg) ^ 0.20 GeV"\ if the b quark mass 
m ^ 4.8 GeV and Tg ^ 0.2 ^ 0.3 fm is taken. Then, if the value of the gluonic condensate 
G2 = 0.05 ± 0.02 GeV^ is used, one finds the estimate in quenched QCD 

cnp < 0.03 ± 0.01 MeV. (11) 

In full QCD an even smaller value is found. This value of cnp is much less than both |aNp|, 
Eq. (7), and the experimental errors. Therefore it can be neglected in the tensor splitting c 
and also in ?]np = |cnp — ^np, i.e. we take here t^np = onp- 

4 Fitting Conditions 

To fit the experimental data 

aexp(lP) = 14.23 ± 0.53 MeV, aexp(2P) = 9.39 ± 0.18 MeV, 

Cexp(lP) = 11.92 ±0.25 MeV, Cexp(2P) = 9.14 ±0.25 MeV, (12) 

the following conditions have to be satisfied: 

atot{nP) = flp^ + flp^ - 7r~^{^'^) = ae^pinP) 

Ctot(nP) = =Cexp(nP). (13) 



4 



As seen from Eqs. (2) and (4), the l.h.s. of these expressions strongly depend on the 
normahzation scale /z, but the combination t] does not. The fitting condition for t] is 

r^(nP) = ^/4 + -^(0=r^e.p, (M) 
irm Im^ 

where the experimental values for r^^xp""* ^'^^ 

Vc:.p{lP) = 3.65 ± 0.9 MeV, ?7exp(2P) = 4.32 ± 0.4 MeV. (15) 
The condition ([I^) does not depend explicitly on and can be rewritten as 

^^f,{nP) = r^exp - = A(nP), (16) 

hence the strong coupling constant can be expressed as 



For a chosen interaction and quark mass m, A{nP) and fi{nP) are known numbers and one 
can immediately determine as(/^)- 

The extraction of as(/i) from the condition (p!?]), in general, extremely simplifies the fit 
and also puts strong restrictions on the possible choice of the normalization scale /x. Just 
this condition was exploited in [|T0| to determine as(/i) for the IP state in charmonium. In 



charmonium rjf^xp ~ 24 MeV and the typical value of A 7 -i- 8 MeV is not small, so the 
uncertainty in the extracted value of ^^(/i) is about 10%. 

In bottomonium the typical values of |aNp| are found to be smaller: |aNp(lP)| = 2.6±0.2 
MeV (see Table D and |aNp(2P)| = 1.95±0.10 MeV (see Table 1). As a result, the numerical 
values of A{nP) to be substituted in Eq. ([T7|) are small: 



A(1P) = 1.05 ± 0.9(exp) ± 0.15(th) MeV, A(2P) = 2.4 ± 0.4(exp) ± O.lO(th) MeV. (18) 

The theoretical uncertainties in this equation are caused by the uncertainty of the value of 
onp in the Thomas interaction. Still, for the 2P state the total error in A(2P) is not large, 
about 20%, and therefore ^^(/i), proportional to v^, can be determined from the condition 
( |17D with an accuracy of about 10%. Our calculations show also that the matrix element 
in Eq. ( P^ is practically constant and therefore the theoretical error in Eq. (|TB|) coming 
from /4 is small. 

For the IP state the experimental error in r^exp, Eq- dH), as well as in A, Eq. (18) 



IS 



large: it comes mostly from the experimental uncertainty in the ^^^(IP) mass. Therefore 
A(1P) can vary in a wide range: < A < 2.0 MeV and the relation (|l^ cannot give an 
accurate value for as(/i). Instead, for the IP state one needs to use the conditions (|T3|) which 
are /i-dependent and less restrictive. 



5 Dependence on Scale 

The second-order perturbative corrections to the spin-orbit and tensor splittings, which are 
not small, explicitly depend on the scale fi. In Eqs. (H, ^) ln(yu/m) enters with the large 
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coefficient 25/6 and therefore the choice /i = m (causing this logarithm to vanish) can give 
rise to inconsistent results. Just this choice was taken in |^ where two 6-quark masses, 
m = 4.6 GeV and m = 5.2 GeV, were analysed. We shall discuss here some results of 
ref. §. 

From the fit in it was obtained that the value ^^(m) extracted from the tensor and the 
spin-orbit splittings are slightly different and for the 2P state this difference is increasing. 
(Here a(/i) or a(/i) denotes the fitted (extracted) value of the strong coupling constant.) 

Also, for the 2P state a(5.2) = 0.26 ± 0.01 is a bit larger than a,(4.6) = 0.25 ± 0.01 
for the smaller 6-quark mass, in contradiction with the standard behaviour of the running 
coupling constant Q;s(g^). 

The extracted value, as{rn) ^ 0.25 -r- 0.27, turned out to be about 20% larger than the 
values as{4:.6) and 0^(5.2) calculated with the conventional value of A^'^^ Eq. ([19|) : as(4.6) = 
0.22 ± 0.01, as(5.2) = 0.21 ± 0.01. 

In the calculations that follow, it will be easy to compare our results with those from 
ref. because in both cases the same perturbative interaction and linear potential ar were 
used. However, the calculations of ref. were done in the nonrelativistic case (for fixed 
a = 0.2 GeV^ and two 6-quark masses). Here both relativistic and nonrelativistic cases will 
be considered and a, m, and aes of the Coulomb potential will be varied in a wide range. 
From our analysis it will be clear that the inconsistencies in the as{f^) behaviour mentioned 
above, are related to the a priori choice n = m made in 0. 

At this point it is worthwhile to note that at present the QCD constant A'-"-'-* is well 
known for Uf = 5, because ^^(M^) = 0.119 ±0.002 is established from different experiments: 
A(5)(2 - loop) = 237li MeV and A'^^\3 - loop) = 219ti MeV are given in 0. Then from 
the matching of as(/i) at the scale fi = rhh {rhh is the running mass in the MS scheme) and 
taking rhb = 4.3 ± 0.2 GeV [0 one can find A('')(3 - loop) = 296t^9 MeV or in the two-loop 
approximation 

A(4) 

IS 

A(^)(2 - loop) = 338i;^? MeV. (19) 

It is of interest to compare Q;s(/i) for A*^^) given by Eq. ([191) with the fitted values as(/i) 
used in different theoretical analyses: 0^(3.25) = 0.251 ±0.009 whereas in [Q the fitted value 
a,(3.25) = 0.33; a,(4.60) = 0.221 ± 0.007 while in § «s(4.6) ^ 0.27. In both fits the 
extracted values appeared to be about 20% larger. 

This 20% difference implies either very large values of A*^^) or an essentially smaller scale 
of yU. For example, as(/io) = 0.33 with the conventional A^^^ Eq. ([T9|) , corresponds to 
/io = l-80±ojg GeV instead of ft = 3.25 GeV in |^ and this fiQ would be in good agreement 
with the one cited in ||^, 0] and with our result (see Section |^). 

In our present analysis when different sets of parameters are taken, we shall impose two 
additional restrictions: 

1) For the given P-state the extracted value of as(/i) must be the same for the tensor and 
the spin-orbit splittings, because both interactions have the same behaviour and they 
also have the same characteristic size (momentum). 

2) Only those sets of parameters for which the fitted two-loop value of (5(/i) corresponds to 
the conventional value of A"^^) in two-loop approximation, Eq. (p!9D, are deemed appropriate. 



6 



6 Static Potential 



In heavy QQ systems the spin- dependent interaction contains the factor and therefore 
it is small and can be considered as a perturbation. For the unperturbed Hamiltonian we 
considered two cases,relativistic and nonrelativistic, 



H§ = 2^p + m^ + V,,{r) (20) 

or 

H^R = ?1 + V,^{r). (21) 
m 

Here a static potential, Vst{r) = V^f (r) + V^^^(r), needs some remarks. The perturbative 
static potential is now known in two-loop approximation |2D], but for our discussion it is 
enough to take it in one- loop approximation from [Q: 

V£ = (22) 

6 r 



Here the vector coupling constant av{r) is expressed through as{n) in the MS scheme in 
the following way [Q: (a<i(/i) << 1) 



l + ^(a. + |(ln(H+7.)^ 



1 - ^ (ai + f (ln(/xr) +7^;)) /5o ln((Aijr)-2) ' 

In Eq. ( |23| ) we have used as{fi) = 47r/[/3o ln(/i^/A|^)], and the conventional QCD constant 
in coordinate space: Ar = A-j^exp^'jE + a) where a = 2ai/ [3q. We see that the dependence 
on 11 disappears. The constants are: (3q = 11 — 2?t,j/3, so for nj = 4, /Sg = 25/3; ai = 
31/12 - 5n//18, so for n/ = 4, ai = 53/36. 

This expression is valid only for small radiative corrections or small distances: re'^'^A^^^ -C 
1 or r ^ 2 GeV~^ = 0.4 fm (A*^^) ^ 0.3 GeV). However, in bottomonium the sizes of the 
different states are varying in a wide range, e.g., typical values of the root-mean-square 



radius, R{nL) = \J {r'^)nL, are: 

R{1S) = 0.2 fm, R{1P) = 0.4 fm, R{2S) = 0.5 fm, R{2P) = 0.65 fm, 

R{3S) = 0.7 fm, R{3P) = 0.85 fm, R{AS) = 0.9 fm. (24) 

These numbers are practically independent of the choice of the static potential parame- 
ters and the confining potential provided the chosen potential reproduces the bottomonium 
spectrum with good accuracy. 

From Eq. ( p^ ) one can see that the sizes of the nL states run from 0.2 fm to 0.9 fm. 
Therefore the perturbative potential, Eq. (p3|), valid for r -C 0.4 fm, can be used only for 
low-lying states. For the IS', 25", and IP states this perturbative interaction (also with 
two-loop corrections) was analyzed in detail in refs. P, and there it was found that (i) 
for the IS* and 25* states the values of fi are different and (ii) is smaller in the 25* state. 



7 



Therefore, one can expect that for every level a specific consideration is needed to determine 
or asin). 

To describe the 2P state, the size of which is about 0.65 fm, or the bb spectrum as a 
whole, a different approach is needed. Here we suggest instead of the perturbative potential 
Eq. (p2D to use the perturbative potential in background vacuum field, VB{r): 



VB(r) = - 



4 aB(r) 



in momentum space 



4 An 
3 ~^ 



2 -'2 

q =q . 



(25) 



(26) 



In this potential ^^(g^) is a vector coupling constant in vacuum background field which was 
introduced in [|T5| and applied to e+e^ —>■ hadrons processes in 
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Pi IntBjq) 
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tsiq) = In 
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m 



B 



A2 



(27) 



with = 25/3. For the vector coupling constant, aviq"^), A differs from A in the MS scheme: 
^« - ASlltl MeV, a = 1 - ^ 



A = A^ e'' 



0.35333, and A^^ was taken from Eq. (|I9D . (In 
the MS scheme Ab and Ajjg coincide for = 4, 5 because of their identical behaviour at 
large [0-) The background mass tub was found from the fit to the charmonium fine 
structure in where = 1.1 GeV was obtained. 

In coordinate space aB{r) can be calculated from the Fourier transform of the potential 
Eq. (H) with aij(g2) given by Eq. Then 



Po 



dq 



sm qr 
qtB{q) 



Pi_ \ntB{q) 
Pi tB{q) 



(2^ 



The strong coupling constant in vacuum background field maintains the property of asymp- 
totic freedom at small r, r ^ A^^ and r ^ ttt^b^, 



0) 



27r 



aB[r 

Here the function 7 = 7(r) is 

7 = = 7£; + S, S 

or at small r 



Po ln(Ae'''r) 



E 

k=l 



-mBr) 



k\k 



(29) 



(30) 



7 = 7s 

whereas in standard perturbative theory 7^ : 
distance r in Eq. (|3lD the expression Eq. (|2^) is always bounded 



rriBr, (31) 
7e = 0.5772. Due to the dependence on the 



For large r , r ^ , the limit of ^^(r) in Eq. (pSf) tends to a constant, denoted as 



0:5(00) and called the freezing value: 

Att 

0^(00) 



Poto 



Pi In to 
Poto 



In^ 

A2 



(32) 
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From the integral Eq. (^) it can be shown that the freezing value is the same in coordinate 
and in momentum space, aB{r ^ oo) = otBi^if = 0). The properties of asir) were discussed 
in 



T0| , |T3|, Q and a detailed analysis of ^^(r) will be published elsewhere. 



In the intermediate region, 0.2 fm < r < 0.9 fm, aB{r) approaches rapidly the value 
q;b(oo). 

Therefore, to study the bottomonium spectrum as a whole it is convenient to introduce 
an effective constant aeff: 

(^sij) = «efr + ^Oisir), cteff = constant, ^ «efr, (33) 

and to consider the contribution from the term 5Vb{t)., 

SVs(r) = -1 (34) 

6 r 

as a perturbation. Then in the Hamiltonian (22) the static interaction 

Voir) = -| ^ (35) 
6 r 

will be taken into account as an unperturbed interaction. 

For the nonperturbative interaction a linear form ar will be taken and therefore the static 
potential in the unperturbed Hamiltonian Vo(r), 

Vo(r) = -^ — + ar + Co (36) 
6 r 

coincides with the well known Cornell potential. Later, the values of the string tension a 
will be varied in the range 0.17 -r- 0.20 GeV^. 

We shall present a detailed analysis of the bb spectrum in a separate paper. 



7 Relativistic Corrections 

There exists the point of view that in bottomonium the relativistic corrections are small 
because of the heavy b quark mass. Indeed, the comparison of levels and mass differences 
for the Schrodinger equation and the Salpeter equation, Eqs. ( ^01 , ^1]) , in general, confirms 
this statement (here the static potential is supposed to be the same in both cases). In Table 
|I| the bb mass differences are given for two typical sets of parameters. 
From Table |I] one can see that 

(i) Relativistic corrections are small for large mass differences like M(n, L) — 
M{n - 1, L) or M{nL) - M{n, L - 1); 

(ii) For close lying levels, like Ai = M{2S) - M{IP) and = M{?,S) - M(2P), 
the corrections are essential, about 15%, and to get Ai and A2 close to the 
experimental data it is necessary to take into account the contribution from 
the perturbation 5VB{r) Eq. (P^. In the relativistic case the influence of the 
phenomenon of asymptotic freedom appears to be more essental than in the 
nonrelativistic (NR) case. 
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Table 1: Bottomonium level differences (MeV) for the Schrodinger and the Salpeter equa- 
tions. 



Mass differences 


Set I, ttes = 0.3545 


Set II, aeflf = 0.36 


Exp. val. 






m = 4.737 GeV 


m = 4.81 GeV 


(MeV) 






a = 0.20 GeV2 {^) 


a = 0.18 GeV2 






Rel. 


Nonrel. 


Rel. 


Nonrel. 




M{2S) 


-M{1S) 


554.34 


551.97 


556.55 


550.03 


562.9 ±0.5 


M(35) 


- M{2S) 


350.43 


354.78 


335.62 


338.49 


332.0 ±0.8 


M(45) 


- M{3S) {^) 


285.93 


291.83 


270.63 


275.30 


224.7 ±4.0 


M(1P) 


-M{1S) 


458.04 


439.66 


473.49 


450.15 


439.8 ±0.9 


M(2P) 


-M(IP) 


359.67 


366.75 


342.55 


348.70 


359.8 ± 1.2 


M(2^) 


-M(IP) 


96.31 


112.31 


83.07 


99.88 


123.1 ± 1.0 


M{3S) 


- M(2P) 


87.06 


100.34 


72.82 


89.67 


95.3 ± 1.0 



This set was taken from ref. p2| 



The 45 level lies above the BB threshold 



Table 2: IP-state matrix elements for the Schrodinger and the Salpeter equations. 



Matrix element 


Set I (") 


Set II (") 




Rel. 


Nonrel. 


Rel. 


Nonrel. 


v/(r2) (GeV-i) 
(r"i) (GeV) 
(r~^ In mr) (GeV^) 
(r-3) (GeV3) 
/4(1P) (GeV^) 


1.994 
0.633 
0.675 
0.551 
0.0876 


2.039 
0.614 
0.631 
0.483 
0.0685 


2.008 
0.631 
0.681 
0.556 
0.0871 


2.054 
0.612 
0.636 
0.485 
0.0673 



For the parameters see Tab. |T] 



Table 3: 2P-state matrix elements for the Schrodinger and the Salpeter equations. 



Matrix element 


Set I C^) 


Set II {") 




Rel. 


Nonrel. 


Rel. 


Nonrel. 


v/(r2) (GeV-i) 
(r-i) (GeV) 
(r^^ In mr) (GeV^) 
(r-3) (GeV3) 
/4(1P) (GeV2) 


3.177 
0.477 
0.495 
0.504 
0.1060 


3.263 
0.455 
0.448 
0.414 
0.0783 


3.235 
0.469 
0.489 
0.496 
0.1025 


3.320 
0.448 
0.443 
0.406 
0.0748 



For the parameters see Tab. |l] 
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The relativistic corrections are becoming essential for some matrix elements, which de- 
termine the fine structure splittings (see Table Q). To calculate them in the relativistic case 
(for the Salpeter equation) the expansion of the wave function in a series over Coulomb-type 



functions was used as it was suggested in . The numbers obtained have a computational 
error < 10~^ (the dimension of the matrices D was varied from D=20 to D=40). 
From the numbers given in Table |^ one can conclude that 

• for IP and 2P states the root-mean-square radii practically coincide in the relativistic 
and the NR cases; 

• for the matrix element {r~^) the difference between both cases is small, about 3% for 
the IP state and about 5% for the 2P state; in the relativistic case and therefore 
\a^p{nP)\ is slightly larger. 

• in the relativistic case the values of (r~^ln mr) are about 7% (10%) larger for the 
1P(2P) state for given set of chosen parameters; 

• for the Salpeter equation the matrix element {r~^) is larger by about 14% (22%) for 
the 1P(2P) state; 

• the largest relativistic correction was found for the factor given in Eq. (^). This 
difference is about 30% for the IP state and 36% for the 2P state. The numbers given 
do practically not change for different sets of parameters. So our averaged value of 
fiijiP) (aeff > 0.35) are: 

/4(1P) = 0.085 ±0.010 GeV^ /4(2P) = 0.106 ± 0.008 GeVl (37) 

The theoretical error in Eq. ( |37D (~ 10%) mostly comes from the variation of the h quark 
mass (in the range 4.6 -r- 5.0 GeV). 

The increasing of f^i^nP) in the relativistic case directly affects the values of ^^(/x) ex- 
tracted from the fine structure data because according to Eq. (|1^) 



TTTD./K \ TiP] 

\ 2U{nPy ^^""^^ = V..AnP) - |aNp(nP)|, (38) 



is proportional to /4 and Q;s(/i) is about 15% smaller in the relativistic case. This result 
obtains both for IP and 2P states. 

Therefore, below we shall use only matrix elements calculated for the Salpeter equation, 
in this way taking into account the relativistic corrections. A last remark concerns the choice 
of the quark pole mass, mpoie = rn which enters the Salpeter equation 0. Here we study 
the spin structure of the Xb mesons determined by the spin-dependent potentials now known 
only in one-loop approximation. Therefore the pole mass of the h quark will be taken also 
in one-loop approximation |[23[| : 



m = mpoie = m{m + \ (39) 
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In Eq. (^) m{m?) is a running quark mass in the MS renormalization scheme, its value 
from III] is m = 4.3 ± 0.2 GeV. Then taking A^"^-* from Eq. ([T9|) one finds m in the range 

4.5 GeV<m<5.0 GeV (40) 

Only values of the mass in this range will be used later in our calculations. 



8 Cis{lj) for the 2P State 

For the 2P state ^^(/i) can be immediately found from the relation (|l^) for the chosen static 
potential with fixed parameters tteS; cr? and m. At first, we shall give an estimate of ^^(/i) 
using the following results: 

1. The nonperturbative spin-orbit splitting |aNp(2P)| depends weakly on the choice of 
the parameters, provided the hb spectrum is described with good accuracy 

|aNp(2P)| = 1.95 ± 0.15 MeV. (41) 

2. In Eq. (plSj) the experimental error of r^exp (2-P) is not large and therefore the difference 
A(2P) Eq. (0) is known with an accuracy of about 20%: 

A(2P) = r/cxp(2P) - |aNp(2P)| = 2.40 ± 0.04(exp) ± 0.15(th) MeV. (42) 

3. In our calculations the matrix element fi{2P) is changing in the narrow range: 

/4(2P) = 0.106 ± 0.008 GeV^. (43) 

Then, from the fitting condition (0) and the numbers given in Eqs. (PD-(^3|) the lower 
and upper bounds of as{f^) can be determined. 




TTl / TTl 

0.37<aJu)<J — 0.46 (44) 
mo V mo 

Here a normalization mass, mo = |M(T(1S')) = 4.73 GeV , was introduced for con- 
venience. Here and below all numbers were calculated in relativistic case, i.e. for the 
Salpeter equation. 



From the estimates ( ^^ it is clear that for the 2P state a;s(yu) ~ 0.40 turns out to be 
large for any set of the parameters of the static interaction. It is essentially larger than 
that obtained in ref. where q;s(3.25) = 0.33 and in ref. |^ where 0^(4.6) = 0.26. In our 
calculations large values of as(/i) are extracted irrespectively to the value of the scale n, 
which is still not fixed. 

However, as{n) in Eq. (^) is varying in a rather wide range and its value is sensitive to 
small variations of the factors entering the condition (p^ . The value of asifJ^) is decreasing 
if the constant aefr of the static interaction is growing. In our numerical calculations the 
values of aes is supposed to be in the range, 

0.35 < aoff < «s(g^ = 0) ^ 0.48 (45) 
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with a b quark mass from the condition (^^. 

With the restriction (|45|) the fitted values of asifi) appeared to lie in the narrower range, 

= 0.40 ± 0.02(th) ± 0.04(exp) (bb). (46) 

Here the experimental error comes from ?7exp, Eq. (p!8|), and the theoretical error is due to 
the variation of aes, "m-, and a. 

In the extracted value as(/i), Eq. (^), the scale /x is still not specified. To find /i2 it is 
better to use the condition c(2P) = Cexp,(P^, for the tensor splitting, because the theoretical 
uncertainty connected with the nonperturbative contribution to c(2P) is negligible, cnp < 
0.05 Mev. This condition (p!3| ) turns out to be satisfied for the scale, 

fx = 1X2 = 1.02 ± 0.02 GeV, (47) 

which has a small theoretical error, 2%, while the extracted value of (5(/i), Eq. (^61), was 
determined with an accuracy of 15%. 

It is of interest to compare a(l.O) ~ 0.40 with the value found in perturbation theory. 
The scale 112 ~ 1.0 GeV is small, less than the running mass of the c quark, = 1.3 ± 0.2 
GeV ||l|], therefore ^^(l.O) should be calculated with A = A'^'^^ (2-loop), Uf = 3. The value 
of A(3) can be found using the matching condition at = rric and the value of A*^^-* (2-loop), 
Eq. (|19D. Then 

A(^)(2 - loop) = 384^1 MeV (48) 
and correspondingly the two-loop strong coupling constant is 

a,(1.0) = 0.53l°|, (49) 

which is 30% larger than our fitted value given by Eq. (|46|). It was suggested in |10[ that 
this decreasing of a{fl) at the scale /i2 = 1 GeV can be explained by the behavior of asifJ^) 
Eq. (^^ in the vacuum background field, thus demonstrating the phenomenon of freezing 
of as{fi) . In |T^, from a fit to the charmonium fine structure, the background mass in 
Eq. ( P7D was found to be (in the MS renormalization scheme) 

TUB = 1.1 GeV, Ag^(2 - loop) = 400l^° MeV (cc). (50) 
Our extracted value of q;(1.0) in Eq. (^) corresponds to the close value of A^, 

Ag)(2- loop) = 4201^° MeV {bb). (51) 
Note also that for the IP state in charmonium the value 

a(l.O) = 0.38 ± 0.03(th) ± 0.04(exp) (52) 
practically coincides with a (1.0) in bottomonium, 

a(l.O) = 0.40±0.02(th) ±0.04(exp), /is = 1.02 ± 0.02 GeV. (53) 

This coincidence is not, in our opinion, accidental: both states, the cc IP state and the bb 
2P state, have the same size: R = \/ < r'^ >np = 0.62 ^ 0.65 fm. This coincidence of the 
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Table 4: Fine-structure parameters for the 2P bb state. 





Set I (") 


Set II (") 


Set III C) 


Exp. val. 


a{fi2) 


0.392 


0.429 


0.386 




fi2 (GeV) 


1.03 


1.02 


1.03 




(MeV) 


-2.62 


-3.14 


-3.35 




ctot (MeV) 


9.12 


9.11 


9.17 


9.1 ±0.2 


onp (MeV) 


-2.12 


-1.83 


-1.80 




a? (MeV) 


17.61 


18.33 


18.77 




a'p^ (MeV) 


-6.12 


-7.19 


-7.52 




atot (MeV) 


9.37 


9.32 


9.45 


9.4 ±0.2 



(") For the parameters see Tab. |l] 

(^) aeff = 0.386, a = 0.185 GeV^, m = 5.0 GeV. 

values of as(/i) and of the sizes indicates that for the bound states the scale /z is characterized 
by the size, but not the momentum, of the system. This result is in agreement with the 
predictions of refs. @, |^. 

With the use of the fitted values as(Ai2), Eq- (0), the theoretical number obtained for the 
spin-orbit splitting atot automatically satisfies the third fitting condition Eq. (pIBf). Calculated 
numbers of a and c are given in Table H for three different sets of parameters. From these 

{2) (2) 

numbers one can see that the second order radiative corrections Op and Cp are negative 
and rather large: about 25% for the tensor and 40% for the spin-orbit splittings. 

Note that we have met here no difficulty to get a precise description of the tensor and 
spin-orbit splittings for the 2P state simultaneously, in contrast to the results of ref. [^, 
where some difficulties have occurred, in our opinion, because of the choice jl = m (see the 
discussion in Section ^. 

9 ci^s(/i) for the IP State 

For the IP state the scale-independent condition (|l^ cannot be used directly, because the 
important factor A(1P) in Eq. (|T7p has a large experimental error. So in this case one needs 
to use the two /i-dependent conditions, Eq. (0), on the splittings a and c. 

There exist a lot of variants where these two conditions can be satisfied. However, in 
many cases the two- loop values Q!s(/ii) and fii, extracted from those fits, correspond to a 
very large value of the QCD constant A'^^^ Therefore, the additional requirement (21) that 
A^'') (2-loop) should have a value in the range 307 MeV < A(4) < 371MeV, is necessary. If 
this restriction is put, then in our calculations the extracted scale /ii appears to lie in the 
narrow range, 

/ii = 1.80 ±0.10 GeV (54) 

and 

a(/ii) = 0.33 ± O.Ol(exp) ± 0.02(th). (55) 
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Table 5: Fine-structure parameters for the IP bb state. 





Set I (") 


Set II («) 


Set III {') 


Exp. val. 


a{fi2) 


0.335 


0.340 


0.32 




/ii (GeV) 


1.80 


1.85 


1.90 




c^P^ (MeV) 


0.96 


1.03 


0.90 




ctot (MeV) 


11.93 


11.92 


11.91 


11.92 ±0.20 


aNP (MeV) 


-2.82 


-2.45 


-2.44 




(MeV) 


16.46 


16.34 


16.52 




a^P^ (MeV) 


0.12 


0.21 


0.05 




atot (MeV) 


13.76 


14.09 


14.12 


14.23 ±0.57 



(") For the parameters see Tab. |l] 

(^) aeff = 0.386, a = 0.185 GeV^, m = 5.0 GeV. 



Our value for the scale /ii turned out to be very close to that determined in refs. , but 
our fitted value of ^^(/ii) is about 15% larger than the one found in ^ where as(3-loop)=0.29 
and A(^)(3 - loop) = 230 MeV (or A^^) (2-loop)= 250l^E] MeV) is smaller than in our fit. 

For the IP state it was also observed that if a large value a = 0.2 GeV^ is taken, then 
it is difficult to reach a consistent description of the tensor and the spin-orbit splittings 
simultaneously. Therefore here, as well as in the charmonium case |jlO|, the values a = 
0.17 -7- 0.185 GeV^ are considered as preferable. Also the choice of a relatively large b quark 
mass, 

mb = 4.75^4.9 GeV, (56) 

gives rise to a better fit. 

The results of our calculations for the IP state are given in Table |5| from which one can 

(2^ (2) 

see that the second order corrections Cp and ap are relatively small, 8% and 15%, but 
still remain important for a fit to the experimental data. Also in all good fits the effective 
Coulumb constant Oeff hes between a(yUi) and as(At2): 

a(/ii) < acs < a(/U2)- (57) 

In our analysis /i2(2P) is less than /ii(lP) and their ratio is almost inversely proportional 
to the ratio of the radii of these states: 



« 1.7 ^1.8; VVhp = 1.6 ^ 1.65. (58) 

This result is in full agreement with the prediction of refs. ^ about the decrease of the 
scale with increasing principal quantum number. 
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10 Conclusion 



The precise experimental data on the masses of Xbj(l-P) and Xbji'^P) give a unique oppor- 
tunity to determine the QCD strong couphng constant at low-energy scales. In our analysis 
of fine structure splittings we found that: 

1. The relativistic corrections which are small for such characteristics as the bb levels, 
radii, and matrix element {r~^), are nevertheless essential for the determination of the 
factor /4(nP), which is inversely proportional to the extracted value of 

2. From a /i-independent analysis of the 2P state, the value as(/U2) ~ 0.40 was extracted. 
The scale /i2 = 1-0 ± 0.02 GeV, determined from the tensor splitting, appeares to be 
practically unchanged for any chosen set of parameters. 

3. The extracted value a(l.O) ~ 0.40 is about 30% lower than the one found in per- 
turbation theory if A^^) = 384;|^ MeV was used. This value agrees with the fitted 
Q!s(1.0,cc) extracted from the analysis of the charmonium fine structure. This result 
can be naturally explained in the framework of background perturbation theory. 

4. The scale fii ~ 1.8 GeV for the IP bb state obtained here agrees with the prediction 
in but corresponds to the larger value A^^^ (2- loop) =338^31 MeV, which gives rise 
to a,(M^) = 0.119 ±0.002. 

5. The preferred values of the pole mass of the b quark are found to be m = 4.7 -h 4.9 
GeV but from the fine structure analysis we could not narrow their range. 

Our results have confirmed the important observation of Yndurain et al. |P, 0] that the 
strong coupling constant is increasing for states with a larger size or larger principal quantum 
number and this fact is essential in many aspects of quarkonium physics. 

One of the authors (A.M.B.) acknowledges a partial financial support through the grants 
RFFI-DFG No.96-02-00088g 
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